Spin-dependent thermoelectric transport coefficients in near-perfect quantum wires 
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Thermoelectric transport coefficients are determined for semiconductor quantum wires witli weak 
thickness fluctuations. Such systems exhibit anomalies in conductance near 1/4 and 3/4 of 2e'^/h 
on the rising edge to the first conductance plateau, explained by singlet and triplet resonances 
of conducting electrons with a single weakly bound electron in the wire (T. Rejec, A. Ramsak, 
and J.H. Jefferson, Phys. Rev. B 62, 12985 (2000)). We extend this work to study the Seebeck 
thermopower coefficient and linear thermal conductance within the framework of the Landauer- 
Biittiker formalism, which also exhibit anomalous structures. These features are generic and robust, 
surviving to temperatures of a few degrees. It is shown quantitatively how at elevated temperatures 
thermal conductance progressively deviates from the Wiedemann-Franz law. 
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I. INTRODUCTION 

One of the main properties of small confined electron 
systems, intensively studied experimentally and theoret- 
ically in the last decade, is the electrical conductance. 
However, other transport coefficients also serve as a sensi- 
tive probe of new phenomena in such systems, such as the 
thermopower of chaotic quantum dots or of atomic size 
metallic contacts Band most recently, anomalies in one- 
dimensional wires [0. Theoretical investigations predict 
in these systems a range of new properties of transport 
coefficients, such as anomalously enhanced thermopower 
in quantum dots due to the Kondo effect Q and, at 
low temperatures, changes in sign together with linear 
thermal conductance violating Wiedemann-Franz law . 
Anomalies in thermoelectric coefficients are also found 
in standard strongly correlated systems: the Anderson 
model 1^], the Hubbard model and the t-J model ||^. 

In this paper, we extend our recent theoretical study of 
conductance anomalies to include thermoelectric effects 
due to a temperature gradient. Anomalies are related to 
weakly bound electron states within the quantum wire. 
In particular, we consider a small fluctuation in thickness 
of the wire in some region giving rise to a weak bulge. 
If this bulge is very weak then only a single electron will 
be bound. We may thus regard this system as an 'open' 
quantum dot in which the bound electron inhibits the 
transport of conduction electrons. Near the conduction 
threshold, there is a 'Coulomb blockade' and we have 
shown that this gives rise to spin-dependent resonances, 
also in an axial magnetic field, for wires of both rectan- 
gular ^ and cylindrical ||l^ cross-section. 

Experimentally, the staircase structure of the con- 
ductance of quantum wires was reported more than a 



decade ago ||Tl[, and more recent systematic investiga- 
tions showed unexpected structure in the rising edge to 
the first conductance plateau Jl2|-|T5t. 

Here we model a quantum wire as in Ref. JlO| and, 
explicitly, we assume a wire of circular symmetry about 
the z-axis with constant potential, V{r, z) = within a 
boundary ro(z) from the symmetry axis and confining 
potential Vq > elsewhere. This geometry is close to 
that of narrow 'v'-groove quantum wires, which also ex- 
hibit anomalies near the conductance threshold jl^]. To 
be definite, we choose parameters appropriate to GaAs 
for the wire and Alj.Gai_2. As for the barrier with x 
such that Vq — 0.4eV, which is close to the crossover 
to indirect gap. The wire width is taken as ro{z) = 
^ao{l + £,cos^ TTz/ai) for \z\ < ^ai and ^0(2;) = ^ag 
otherwise, i.e., a wire of width oq with a single bulge 
of length ai and width (1 -|- f)ao, as shown in insets of 
Fig. 1(c) and Fig. 2(c). 



II. CONDUCTANCE 

We consider the interacting electron problem with the 
above wire thickness variation in a range which ensures 
that only one electron occupies a bound state and that 
restriction to a single channel near the conduction edge 
is an excellent approximation. This is always the case for 
a very weak smooth variation, i.e. a near perfect wire. 
From numerically exact solutions of the two-electron 
scattering problem, the conductance is calculated from 
our generalisation of the usual Landauer-Biittiker (LB) 
formula ||l^ , to include spin-dependent scattering ||l^ of 
conduction electrons from the single electron bound in 
the potential well. This gives G'(/i) = Go where 



1 



Go = 2e^//i, is the Fermi energy and the transmitivity 
is a weighted average over singlet and triplet channels 

"Ml, 



where 



T{^,) = \%{^,) + \%{^i). 



(1) 



At elevated temperatures we use the LB finite tempera- 
ture extension 



G(/i) = Go 



de 



T{e)de, 



(2) 



where f{e,fi,T) = (1 + exp[(e — fi)/kBT])'^ is the usual 
Fermi function which describes the thermal distribution 
of electrons in the leads. G(/i) is shown in Fig. 1(a) 
and Fig. 2(a) for a wire with relatively small and a larger 
bulge, respectively. Here the energy is measured from the 
threshold of the conductance. As discussed in Ref. pO[ , 
the weak bulge in the wire is equivalent to a shallow po- 
tential well in a perfectly straight wire and if the length of 
the bulge region is small, this effective potential well can 
only accommodate one bound state with the consequence 
that only a singlet resonance in G exists, as observed, for 
example, in Ref. Conversely, if the bulge region 

is longer, both, singlet and triplet resonances contribute. 
For even longer bulge regions with a very shallow effective 
potential well (near perfect wire) , the singlet resonance is 
pushed to lower energy and therefore becomes extremely 
narrow. In this regime, only the broader triplet can be 
resolved at finite temperature |9|, [l0|, a s observed experi- 
mentally in clean gated structures [|l2|-^ . 



III. THERMOELECTRIC EFFECTS 

The LB approach can be extended to include elec- 
trical and heat currents through a region between two 
leads with different temperatures and chemical poten- 
tials |l9|j2^. With T + AT, fi + eU for the left lead and 
T, /i for the right lead, we get 



JQ 



j^^J Af{e)Tie)de, 

(e-M)A/(e)r(e)de, 



and 



A/(e) = /(e, ^l + eU,T + AT) - /(e, /i, T). 



(3) 
(4) 

(5) 



In the linear response regime of vanishing AT and U the 
currents simplify to 
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KM 
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FIG. 1. (a) Electrical conductance G(/i), (b) thermopower 
S(/i), and (c) thermal conductance k(/i) for wire parameters 
ao — 10 nm, ai = 30 nm, ^ = 0.18 and screening length 
p — 100 nm. Other parameters and the numerical method 
is as in Ref. 10. The dashed line in (c) represents Wiede- 
mann- Franz law result for T — 4K. The traces for different T 
are offset vertically for clarity. 



A. Thermopower 

The Seebeck thermopower coefhcient S measures the 
voltage difference needed to neutralize the current due 
to the temperature difference between the leads. In the 
linear response regime the thermopower is given by. 



> 

0) 



03 



CM pa 




5(m)- — ---^ 



(9) 



as is for various systems discussed in Refs. [^^1). Eq. ||is 
formally the same as the Mott- Jones formula for simple 
metals pl| and generalized for a system with stronger 
electron-phonon interactions in Refs. |22]. 
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FIG. 2. As Fig.l, but for longer bulge with parameters 
ao = 10 nm, ai — 60 nm and ^ — 0.1. 

In Fig. 1(b) the thermopower of a narrow wire with 
a small bulge is presented for the same range of tem- 
peratures as G{/j,). Such a result is expected, e.g., for 
the system studied in Ref. ||l5|. The structure reflects 
the singlet resonance observed in the conductance and is 
smeared out at temperatures comparable with the width 
of the resonance. In a wire with small thickness variation, 
but with a longer bulge, triplet resonance scattering also 
exists, as shown in Fig. 2. In the thermopower curve of 
Fig. 2(b), the dominant structure at lower temperatures 
comes from the singlet resonance, though the triplet res- 
onance is still clearly discernible. At higher temperatures 
the triplet structure is washed out first, in contrast to the 
conductance result. Fig. 2(a). At low temperatures only 
the transmitivity at energies close to the chemical po- 
tential contributes to the above integrals and the general 
result Eq. ^ can be related to the temperature dependent 
G{fi) by the following expansion 



TT^kB T (dlnGiji) 



3e 



9/1 



15 G{^l) dfi^ J ' 



(10) 



Our results were calculated using the exact relation Eq. ||. 
However, the leading term in Eq. is a reasonable 
approximation for energies above the singlet resonance 



and up to temperatures where the structure is thermally 
smeared out. This is shown in Fig. 3(a) where we present 
a comparison of S{ii) for the exact result with the ap- 
proximations to first and second order. We see that at 
energies below the resonance, both the linear and cu- 
bic approximations deviate significantly from the exact 
result, Eq. ^. In this regime the conductance is itself 
very small and hence G{fi)~^ d^G{fi) /dfi"' is prone to er- 
ror making calculations and experimental data analysis 
based on this expansion unreliable. 

The thermopower of one-dimensional wires has been 
measured 24 1 and more recently, further anomalies 
related to '0.7 anomaly' in conductance were reported 
The authors of Ref. observe a dip in S'(yu) at en- 
ergies corresponding to the anomaly in G(/i). However, 
the logarithmic derivative with respect to the gate volt- 
age of the measured G exhibits a much deeper minimum 
than the dip in the measured S{ii), which remains well 
above zero even at the lowest temperatures . This clearly 
shows that a simple non-interacting formula is not valid 
in this low temperature regime. Apart from the small 
corrections to the logarithmic approximation to S, our 
model and its solution within the LB framework are 
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FIG. 3. (a) Thermopower as obtained with Eq. ^ for 
T = 2K and parameters used for Fig. 1 (full line) . Dashed and 
dotted lines correspond respectively to the result of Eq. |l^ 
and the linear T approximation (first term in Eq. (b) 
Thermal conductance - parameters as in (a). 

in agreement with the findings of Ref. [^. That is, the 
calculated thermopower is in good agreement with exper- 
iment except at low temperatures where we also predict 
a deep minimum. This discrepancy at low-temperatures 
may well be a many-body Kondo-like effect contained 
within our model but not within the two-electron ap- 
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proximation we have used here and in our earher pa- 
pers. We expect the two-electron approximation to break 
down at low temperatures for which the underlying ex- 
tended Hubbard model, which is the starting point of 
our approach, can be mapped onto a generalised Ander- 
son model with coupling terms that are strongly energy 
dependent | |25||. The standard results for the single impu- 
rity problem ||6[ cannot be applied directly to this effec- 
tive model, which is the subject of current research [ p6[ . 
At very low temperatures, a Kondo-like resonance is ex- 
pected [^, for which many-body effects would dominate 
with a breakdown of formula Eq. ||. 



B. Thermal conductance 

The linear thermal conductance is the heat current di- 
vided by the temperature difference between the leads 
when the chemical potentials are adjusted to give no elec- 
trical current. From Eqs. we see that this is related 
to r(e) by. 



(11) 



For low temperatures this simplifies to Wiedemann- Franz 
law, first term in 



3e2 
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In Fig. 1(c) and Fig. 2(c) k(/i) is shown for T from 
0.2K to 4K, calculated from Eq. |ll|. Comparison of 
Figs. 1(a), 2(a) with Figs. 1(c), 2(c) shows good agreement 
with the Wiedemann-Franz law at lower temperatures 
but there is increasing deviation at higher temperatures 
in the resonance region. For comparison, the dashed lines 
in Fig. 1(c), Fig. 2(c) show the corresponding linear ap- 
proximation result, Eq. |l^. This is also seen in the plot 
of K for T — 2K shown in Fig. 3(b). One of the most 
striking features of these plots is that K(/i), calculated 
from Eq. ^l|, exhibits an anomaly at higher energies than 
the corresponding anomaly in conductance, a prediction 
which is open to experimental verification. 



and G ~ IGq is a spin effect, being a direct consequence 
of the singlet and triplet nature of the resonances. The 
probability ratio 1:3 for singlet and triplet scattering fol- 
lows directly from this and as such is a universal effect, 
not only for conductance but all thermoelectric transport 
coefficients. A comprehensive numerical investigation of 
open quantum dots using a wide range of parameters 
shows that singlet resonances are always at lower ener- 
gies than the triplets, in accordance with the Lieb-Mattis 
theorem for bound states |^ . 

Thermopower plots show anomalies, related ultimately 
to the Coulomb interaction between a localised electron 
and the remaining conduction electrons. We have shown 
that the lower-energy singlet anomalies in thermopower 
are more pronounced. These should be clearly observ- 
able in wires which show the corresponding conductance 
anomalies, such as the narrow 'hard confined' wires re- 
ported in Ref . |l5|] , or in gated quantum wires under high 
source-drain bias where the singlet anomaly is clearly ob- 
served p^ . 

Finally we conclude by emphasising that although our 
model of a quantum wire with a weak bulge may appear 
rather specialised, it is actually quite general since the 
weak bulge is mathematically equivalent to a weak po- 
tential well in an otherwise perfect wire. As with our 
previous work, we have not investigated in detail the ac- 
tual causes of such weak effective (or real) potential wells 
but point out that they may well be due to quite different 
sources in different experiments, e.g. thickness fluctua- 
tions, remote impurities or gates, electronic polarisation, 
or some other more subtle electron interaction effect. The 
main point is that because the effective potential well is 
shallow, it will bind one and only one electron. The uni- 
versal anomalies in conductance and thermopower are a 
direct consequence of this and occur for a wide range of 
circumstances in almost perfect quantum wires. 
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IV. SUMMARY 



In summary we have, within the framework of the LB 
approach, calculated thermal transport coefficients for 
near-perfect quantum semiconductor quantum wires, ex- 
tending our earlier work on spin-dependent conduction 
anomalies. These anomalies are a universal effect in one- 
dimensional systems with very weak longitudinal confine- 
ment. The emergence of a specific structure G(/i) jGq 
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